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We present a computational and theoretical study of a many-body Brownian ratchet, in which 
a "gel" of multiple, stiff polymerizing filaments pushes a diffusing obstacle. Our results show that 
steady-state dynamics of this system are strongly influenced by a layer of depleted fllament density at 
the obstacle-gel interface. Inter-filament correlations within this molecule-thick layer have dramatic 
consequences for the velocity and structure of the growing gel. These emergent behaviors can be 
captured by mean field theories that emphasize the non-additivity of polymerization forces and 
indicate a key role for the fluctuating gap between gel and obstacle. 



Living systems have evolved many processes that ex- 
ploit fluctuations at the sub-cellular scale to transmute 
chemical energy into mechanical work. These processes, 
generally referred to as Brownian ratchets, propel cell 
motions such as crawling, phagocytosis, cytoplasmic 
streaming, organelle transport, and chromosome separa- 
tion during anaphase |lL| . Their essential physics is illus- 
trated by a two-body model of a polymerization-based 
ratchet proposed by Peskin, Odell, and Oster [5^. Here, 
a single, stiff chain-molecule, such as an actin filament, 
pushes an obstacle with diffusion constant D, such as the 
cell membrane, microbread [iT, or atomic force micro- 
scope cantilever [4j, by polymerizing against it. When- 
ever a spontaneous fluctuation in the obstacle's position 
opens a sufficiently large gap, a monomer of size a can 
stochastically intercalate and polymerize onto the flla- 
ment at an average rate fco„. Because the fllament is 
imagined to be anchored within a rigid gel, and because 
polymerization is taken to be irreversible, this stochastic 
event ratchets the obstacle's diffusive advance, inducing 
a mean obstacle drift velocity v. The master equation for 
this siTi^/e-fllament model is exactly solvable for steady 
states at all values of D and external load / applied to 
the obstacle. In the limit D ^ fco„a^ the resulting f-v 
relation has a simple form 

w afco„cxp(-/a/A;BT), (1) 

where k^T is the available thermal energy [2]. 

The more biologically relevant problem of a many- 
fllament ratchet is considerably less tractable, in part 
because inter-fllament correlations make analytical treat- 
ment more difficult and in part because long relaxation 
times make computational study taxing. The single- 
fllament model can be viewed as a caricature of such an 
A^-fllament system in which different fllaments are as- 
sumed to polymerize and push independently, each bear- 
ing an effective load f /N j5j [6]. According to this ap- 
proximation we could estimate the steady-state drift ve- 
locity V for an unloaded 7V-fllament ratchet by setting 



/ — fdi-ng/N — vk-gT/ DN, the portion of the drag force 
borne by each fllament, in Eq. (111. We then obtain 



in the limit DN ^1. In Eq. ([2| and below we adopt 
units of length and time such that a — 1 and fcon = 1- 

The assumption that the pushing forces of fllaments 
simply add has been used as a basis for many-body sim- 
ulations [7 , and to extract predictions for TV based on 
/ and V from phenomenological models in which N does 
not directly appear (S) E] • It is also inherent to contin- 
uum models of growing gels |10|, lll| wherein the poly- 
merization force on an object results from the integral of 
a polymerization stress governed by a local constitutive 
law. All these approaches in effect assert that dynam- 
ics of the many-fllament problem is well represented by 
growth of a single fllament under an evenly distributed 
load; however, this assertion has never been verifled for 
a many-fllament model incorporating obstacle diffusion 
and stochastic polymerization with microscopic detail. 

In this Letter we report on a set of simulation re- 
sults and approximate theories of the iV-fllament ratchet 
that overcome previous difficulties, allowing scrutiny of 
the routinely neglected correlations among different flla- 
ments' growth. We show that these correlations can in 
fact shape the steady-state velocity and internal struc- 
ture of a many-body ratchet. In particular, our results 
highlight the signiflcance of a microscopic layer of de- 
pleted fllament tip density at the growing edge of the gel. 
Because of this depletion, apparently weak inter-fllament 
correlations are amplifled in their effect on structure and 
growth kinetics. This depleted layer has been discussed 
for related models [TJ [T5J [T3], but its signiflcance, na- 
ture, physical and mathematical origins, and generality 
for A^-body ratchets were not explained. 

The model we study, sketched in Fig. Q, is a gener- 
alization of the A^ = 1 ratchet [51 to the multi-fllament 
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Figure 1: Monomers with height a = 1 (gray rectangles) 
stochastically polymerize at a mean rate kon = 1 onto existing 
filaments with heights Xi, while the obstacle (thick black line) 
at y executes unbiased, one dimensional diffusion. Filaments 
are straight, parallel, rigid, and randomly aligned. The coor- 
dinate s and distances Zi are measured relative to y. Coordi- 
nates u and w are measured relative to X = maxi=i_..._jv Xi. 



case. Like that model ours lacks the chemical intrica- 
cies of any specific biological ratchet, instead focusing on 
generic physical features that are essential to its function. 
In detail, we consider N parallel, straight, rigid filaments 
(the "gel") that push, via polymerization, against a dif- 
fusing obstacle with diffusion constant D and no external 
load. The position of filament tip i advances stochasti- 
cally in the +x direction as a conditional Poisson process, 
taking discrete steps of size a with mean rate kon so long 
as monomer addition would not penetrate the obstacle 
at position y. The obstacle in turn diffuses freely with a 
reflecting boundary condition at the leading edge of the 
gel, X = maxi Xi . We imagine that the base of each fila- 
ment is firmly anchored (as in the N = I model), so that 
Xi increases in time on a fixed one-dimensional lattice. 
Because the actin gels we have in mind are highly disor- 
dered materials, we take the offsets among these lattices 
to be randomly distributed. The number of filaments N 
and the diffusion constant D are the only dimensionless 
parameters in our model. Biologically relevant values of 
D range from D w 0.01 for L. monocytogenese in cy- 
toplasm [5j to Z? « 10 for patches of membrane |14| . 
We consider values of N between = 1 and N = 2^^ 
(= 16,384). 

We generated stochastic trajectories of our model using 
the continuous time Monte-Carlo method, which sam- 
ples dynamics efficiently and exactly up to numerical 
round-off error. In our application of this method, time 
advances stochastically from the most recent attempted 
polymerization event at t to the next attempted event at 
t + T. The random waiting time r between polymeriza- 
tion attempts is distributed as p(r) = 7Vexp(— A^r). At 
each attempt, a randomly selected filament polymerizes 
provided Xi{t) -|- 1 < y{t ~\- r). Between polymerization 
attempts, the obstacle moves diffusively with a reffect- 
ing boundary condition at X{t). We initialized simula- 
tions by setting y — a and randomly assigning 0:^(0) in 
the interval [0, a], and then computed ensemble statistics 
by sampling configurations at uniform time intervals fol- 
lowing establishment of a steady state. Approach to the 
steady state, however, becomes very slow as v approaches 




Figure 2: Dimensionless velocity, v, in steady state as a 
function of the number A'' of pushing filaments. Results are 
shown for simulations (points (•) with thick solid lines to 
guide the eye) and several predicitons: based on additive poly- 
merization forces (dash dotted line, by solving Eq. ([T| with 
/ = v/DN), from FMFT [Eq. 0] (dashed lines), ar id from 
GMFT assuming a = [?(1)/L>]^'^> 1 [Eqs. ([sj and ( |l2|l4| ] 
(thin solid lines). Colors indicate different obstacle mobili- 
ties: L» = 10 (black), 1 (magenta), 0.1 (red), and 0.01 (blue). 
(Inset) 1 — v{N) (colors and markers as in the main panel) 
compared with 1/A'' (solid black line). 

one. To verify that our statistics are representative of 
terminal steady states, we confirmed that several exact 
conditions [detailed below in Eq. ([4])[ were satisfied. 

Our key numerical results include v{N,D) and the 
spatially resolved average density of filament tips. The 
small- and large- A^ limits of v{N,D) are dictated, re- 
spectively, by the Peskin-Odell-Oster solution for A^ = 1 
and the average unobstructed polymerization velocity 
V — 1. Our simulation results presented in Fig. ^ 
show that the crossover between these limits is gradual, 
with V « log A'^ over a large, intermediate range of A^. 
For AT > 1 we find that w sa 1 - S,iD)/N, where i{D) is 
a dimensionless function of D. This asymptotic scaling 
with A^ is consistent with Eq. ([2|, though the prediction 
that £,{D) = D^^ is incorrect. Except for the largest 
value of D we consider, the number of pushing filaments 
required to sustain a given velocity is significantly under- 
estimated by this heuristic prediction (by several orders 
of magnitude when D < 1 and N ^ 1). 

As we show in Fig. ([3^), for A^ > 1 the linear density 
of filament tips at distance s from the obstacle, p{s), is 
marked by a monomer-thick layer of depleted filament 
density adjacent to the obstacle. Within mean field the- 
ories discussed below the many-body feature underlying 
this non-monotonic density profile is the steadiness of 
induced drift when multiple filaments are pushing. Poly- 
merization stalls in s < 1, which for A^ = 1 transiently 
nullifies the effective force of ratcheting. For A^ > 1 stall 
of any one filament does not disrupt collective polymer- 
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ization forces, so that the flux of stalling filaments into 
s < 1 is balanced by obstacle drift (rather than by ob- 
stacle diffusion as for iV = 1). This balance produces a 
microscopic zone of depleted tip density, which amplifies 
the influence of weak inter-fllament correlations. 

To clarify the relationship between steady state kinet- 
ics and microscopic structure at the gel's leading edge, we 
develop approximate analytical solutions to the master 
equation for the time-dependent configurational proba- 
bility P{y, xi, . . . , xn] t) of our iV-filament model |15| . 
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[coordinate-name subscripts (i.e. t, y, s, and u) de- 
note derivatives]. The first term on the right hand 
side of Eq. ([s]) represents free diffusion of the obstacle. 
The remaining terms, which involve the shift operator 
Tif = f{y, . . . ,Xi — 1, . . . , XN]t) and the Heaviside func- 
tion 9(x), represent stochastic growth of the filaments 
as constrained by the obstacle. Because obstacle diffu- 
sion is a continuous process, the mutual impenetrability 
of the obstacle and gel requires the boundary condition 
Pyi{y — ^i}) = to prevent flux of probability to con- 
flgurations that violate constraints of volume exclusion. 
For = 1, Eq. (|3| is identical to the equation studied 
in Ref. H,. 

Steady-state solutions to Eq. ([3| are most naturally 
described in a moving frame. Any coordinate that trans- 
lates with average velocity v can serve as the origin of 
this frame. We will make use of two different choices, 
namely, the obstacle position y and lead fllament posi- 
tion X. While the choice is physically immaterial, it 
determines which kinds of correlations can be treated 
straightforwardly when > 1. Focusing flrst on the 
obstacle's motion, we deflne a set of relative coordinates 
Zi = y ~ Xi. Because the lattice offsets are uniformly dis- 
tributed, dependence of P on the Zi variables should be 
insensitive to the value of y at steady state. 

We derive several exact relationships by evaluating mo- 
ments of Eq. ([3]). In particular, at steady state the 
average y-current (J yPtdx^ dy) and average Xi-current 
(J XiPtdx^ dy) both yield the mean drift velocity v. In- 
tegrating by parts, we can write these two moments in 
terms of average structural properties, speciflcally the 
filament tip density p{s) = /{^.^q} X^i '^(^ ~ Zi)Pdz^ 
at a distance s from the obstacle, and the distribution 
tp{u) = J 6{u — y + X)Pdx^dy of gap distance m = — 



v = Dp{0) = Dm = j^ 



p{s)ds. 



(4) 



The last expression for v in Eqs. Q was used in [2_ 
for the = 1 case; we primarily make use of the other 
two. According to the first relation, v = Dp{0), the effec- 
tive force exerted on the obstacle, v/D, is proportional 



to the average number of fllaments in contact with it. 
The mechanical reasoning behind Eq. ([2| does not ac- 
knowledge this sensitivity of polymerization forces to the 
gel's structure at contact, which is strongly shaped by 
multi-fllament correlations. 

We derive an equation for p{s) from Eq. ^ by multi- 
plying both sides by the density operator 6{s — y + Xi) 
and integrating over xi, . . . , xn , 

Dpss - Dpf'^ (0, s) + Tip - pQ{s - 1) = (5) 

The corresponding boundary condition, Ps(0) — 
^'■^■'(0,0), can be obtained in similar fashion. These re- 
sults involve, but do not determine, the two-point corre- 
lation function p^'^\s, s') = /{^^>o} '^i^j ^ Zi)5{s' - 

(Sj) we introduce a mean fleld 
ects correlated fluctuations in 



Zj)Pdz^. To close Eq 
approximation, which neg 
density at different positions s and s'. 



N -1 
N 



p{s')p{s) 



(6) 



|17| . This approximation, which we refer to as fllament 
mean fleld theory (FMFT) and indicate with an over- 
bar, is conceptually and mathematically distinct from the 
assumption of additive polymerization forces underlying 
Eq. Q. Substituting Eq. (|6| into Eq. ^ yields a 
self-consistent equation 

Dp,, - DTps is) + p(s + 1) - pe(s - 1) = 0, (7) 

where T = p{0). The corresponding boundary con- 
dition is ps{0) = Fp{Qi). The solution of Eq. ([7| is a 
piecewise linear combination of exponentials in s. 

The nonlinear mean-fleld equation (|7| is an effective 
one-body model for the many-body ratchet, which treats 
explicitly the stochastic growth of a single tagged flla- 
ment and replaces the remaining A'^ — 1 fllaments with 
a constant effective driving force J- . Its form is identi- 
cal to the steady-state A^ = 1 master equation ([3]) for 
a system with an external load P , but with p replacing 
P and with the alternative normalization pds — N . 
FMFT determines the mean-fleld driving force F self- 
consistently, requiring that it match the polymerization 
force v/D emerging from the effective one-body problem 
on a per-fllament basis, ? /{N — 1) = v/DN. 

Predictions of FMFT for v and p, obtained numerically 
and shown in Figs. (2|3), improve upon the prediction 
heuristically derived from Eq. ([T| over the entire range 
of A^ and D we have explored (greatly so for large A^ and 
small D). In particular, FMFT captures qualitatively 
the depletion of fllament tip density at the gel's leading 
edge. 

We can estimate p{s) asymptotically when D ^ 1 and 
A^ ^ D^^ (so that DT ^ v ^ 1, see below), shedding 
light on the depleted layer's origin. For s > 1, because 
obstacle drift nearly keeps pace with free polymerization. 
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Figure 3: Steady state gel structure for D = 0.1 and A'^ = 
600, as determined from simulation and mean field theories. 

(a) Filament tip density p(s), from simulation (red), from 
FMFT (dashed green Une), and from GMFT (solid green line). 

(b) p(s) (red line) and p*(s) [Eq. ([8|] (blue line), both from 
simulation, (c) Gap distribution ijAu) from simulation (red) 
and from GMFT (green) (d) p*(s) — p(s) from simulation 
(red) and from GMFT (green). (a,c,d) For GMFT we assume 
a » 1. 



p varies slowly on the monomer scale a = 1. Using the 
gradient expansion p{s + 1) « p[s) + ps{s) + \pss{s), we 
find that the mean drift of filament tips, (1 — DT)ps, 
is balanced by fiuctuations in polymerization, ^Pss, ob- 
taining p(s) p{l) exp[— 2(1 — v){s — 1)] for s > 1. From 
this form and the third of Eqs. ^ we also compute 
£,{D) ~ 5P(1)- [In terms of p{l) and v, these predic- 
tions of FMFT agree very well with simulation.] For 
s < 1, polymerization fiuctuations are quenched, and the 
dominant balance of Eq. ^ is instead between obsta- 
cle drift, DJ^ps, and the approximately constant infiux 
of stalling filaments, p{s -I- 1) ~ This balance, the 

boundary condition for Eq. ([7|, and the first of Eqs. ^ 
immediately yield p{l) ~ D^p{Q) ~ D~^, and predict 
p{s) ^ D^^ + sD^^ for s < 1. The actual solution to 
FMFT is smooth; here we neglect an exponential transi- 
tion layer near s = 1 with width D. The general shape 
of p from FMFT agrees with simulation, but v and p(l) 
differ significantly from simulation when N ^ Z?~^ ^ 1. 

We checked the factorization hypothesis of FMFT, Eq. 
^ , using simulation data. Within FMFT the quantity 



~ - 1 p{s) 



(8) 



should equal p{s). and is indeed barely distinguishable 
from p{s) except within the stall region s < 1 [see 
Fig.([3|3)[. Steady-state growth, however, is determined 
by the filament density precisely at contact, which is 
much lower than p{l) for iV ^ 1 due to the depleted 
layer. Relative to /o(0), deviations from Eq. ([6| become 
in fact significant as s — >■ 0, leading to significant discrep- 
ancies between v and simulation results. 

One source of correlations neglected by FMFT is the 



transient, continuous motion of the obstacle between 
discrete polymerization events. This motion uniformly 
shifts the instantaneous filament density profile as a func- 
tion of the relative coordinate s. The resulting change in 
density at different values of s is highly correlated, involv- 
ing all N filaments, and is difficult to capture in a theory 
based on the co-moving frame of the obstacle. Viewed 
relative to the lead filament, however, the gel's instanta- 
neous density profile is insensitive to transient obstacle 
diffusion. This simple fact motivates a change of refer- 
ence frame, in which the fluctuating variables include the 
gap u between gel and obstacle, and the lag Wi = X — Xi 
of filament tips behind the lead filament. 

Parallel to the development of FMFT, we compute 
from Eq. ([s]) an exact equation for the gap probability 
distribution tpiu), 



D^uu = / [n( 

^0 



u, w)Q{w + M — 1) — n(M — w + 1, w)] dw, 

(9) 

together with the boundary condition V'm(O) = 0. Solving 
Eq. ([9]) for i^{u) requires knowledge of the joint probabil- 
ity n(u,w) EE jY.'^^i5{u-y + X)6{w-X + Xj)Pdx^dy 
characterizing correlations between gap size u and fil- 
ament density fiuctuations at a distance w behind the 
lead filament. As before, we construct a closed set of 
equations through an approximate factorization. 



(10) 



here asserting statistical independence between fiuctua- 
tions in gap size and those of the gel's internal struc- 
ture. This alternative mean field hypothesis, which we 
call gap mean field theory (GMFT) and indicate with 
tildes, focuses attention on a distinct description of aver- 
age gel structure, (j{w) = J J2^=i^{'^ ^ ^ + d.y, 
in which filament density is spatially resolved relative to 
the lead filament position. 

Substituting Eq. (10 1 into Eq. ^ yields a self- 
consistent equation 

Dipuu — / ^/j{u)Q{w + u — 1) ip{u — w + 1) (j{w)dw 
Jo L J 

(11) 

that requires as input the average filament density a{w). 
In our approach this input is obtained by solving the 
equations 



Pis) 



ip{s — w)(f{'w)dw 



p(2)(0,s) = ij{0) lim a{s + e) 



e-yO+ 



(12) 
(13) 



together with Eq. ([s]). Eqs. ( T2]|T3 ) follow from the exact 
relations 

/•OO 

/3(s) = / Il{s — w,w)d'w 
Jo 

p(^'(0,s) = lim n(0,s) 

e->-0+ 
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using the factorization hypothesis Eq. (10 1. 

GMFT, given by Eqs. and pl|T3l ), simplifies when 
the polymerization of filaments into the gap is fast com- 
pared to the diffusion of the obstacle over the distance 
a = 1. In this limit the average gap size, i.e., the width of 
the distribution ip{u), is much smaller than a; the scale 
of variation of a{w), on the other hand, cannot be signif- 
icantly smaller than a, due to the random offsets among 
different filaments' growth lattices. Treating a{w) as con- 
stant over the scale of decay of tp{u), we can reduce Eq. 



(Ill to 



'ip{u')du' , 



(14) 



where ii = au, a = [(t{1)/ D]^^"^ ^ 1, and — 
^{u/a)/a. For u 3> 1, i^{u) ^ cexp{—u~^^^) for a con- 
stant c, in contrast to the exponential distribution of gap 
sizes for A'^ = 1. The correspondingly more rapid de- 
crease of ^ with increasing u amounts to a renormaliza- 
tion of the diffusion constant D in the context of FMFT. 

To determine the renormalized valu e o f D, we first 
computed ip{u) numerically from Eq. (14 1. Use of this 
result in Eqs. ( 12p3 1, Laplace transforms, and the sep- 



aration in the scales of variation for -0 and a yields the 
gradient expansion 



p(0)p(s) + xPs 



(15) 



where x = ij{0) u'ip{u')du' w 0.6844. The first 
term in Eq. ^ differs from Eq. ^ of FMFT only 
in the absence of the pre-factor {N — l)/N. The sec- 
ond term in Eq. (15 1 is an 0(1) renormalization of 



D due to correlated density fluctuations neglected by 
FMFT. This renormalization alters the boundary con- 
dition for p at s — 0, resulting in significant improve- 
ments to p and V as compared to FMFT. Asymptoti- 
cally for N > £1-2 > 1^ GMFT predicts p(0) - D-^ 
and psiO) /5(1), as in FMFT. It adjusts p(l), how- 
ever, p{l) - (1 - x)~^D-^ « 3.17p(l) and yields |(D) = 
i(l - xy^D-^ ~ 3.17^(0). We also find a{l) - p(l), 
and so a (1 - x^^^^D^^ ~ lAlD^^ as iV ^ oo when 
D ^ 1. In Fig. (|3|i), we compare Eq. (151 with sim- 
ulation. The solution to Eq. (14 1 agrees very well with 
simulations even for a ~ 10 [see Fig. (|3|:)]. Eq. (13 1 also 
agrees well with simulation. 

Note that like FMFT, GMFT is exact for TV = 1. Of 
the two mean field theories, GMFT agrees more closely 
with simulations when iV 3> 1. As such, GMFT is the 
natural starting point to analyze other physical effects in 
A^-filament ratchet systems such the application of a load 
to the obstacle, filament and obstacle elasticity, Arp2/3 



mediated filament-branching, three-dimensional gel "ar- 
chitecture", and adhesion between the obstacle and gel. 
Our results here suggest that many-filament ratchets are 
especially sensitive to physical effects that modify the 
dynamics of the fluctuating gap, but none of the afore- 
mentioned physical effects alter the applicability of the 
factorization hypothesis underlying GMFT. By separat- 
ing the fast dynamics of the gap from the much slower 
chemical kinetics and elastic fluctuations of the gel, the 
structure of GMFT both correctly handles many-body 
correlations in growing gels and naturally lends itself to 
multi-scale simulation, which can overcome the problem 
of slow relaxation in more complex systems. 
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